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A central aim of physics is to describe the dynamics of physical systems. Schrodinger's equation does this for 
isolated quantum systems. Describing the time evolution of a quantum system that interacts with its 
environment, in its most general form, has proved to be difficult because the dynamics is dependent on the 
state of the environment and the correlations with it. For discrete processes, such as quantum gates or 
chemical reactions, quantum process tomography provides the complete description of the dynamics, 
provided that the initial states of the system and the environment are independent of each other. However, 
many physical systems are correlated with the environment at the beginning of the experiment. Here, we 
give a prescription of quantum process tomography that yields the complete description of the dynamics of 
the system even when the initial correlations are present. Surprisingly, our method also gives quantitative 
expressions for the initial correlation. 

There is a rich history to the studies of decoherence of quantum systems due to the interactions with the 
surrounding degrees of freedom. When the dynamics of the system (<S) is Markovian it can be described by a 
master equation 1 " 3 . Nowadays many researchers are interested in systems that are non-Markovian, as there 
is mounting evidence that some natural systems of importance may be non-Markovian 4 and such features may 
allow to manipulate and control quantum systems in desired ways. There is also a great deal of interest in sy- 
stems that are initially correlated with their environments (£) because non-Markovianity and initial system- 
environment (SS) correlations are intimately related 5 " 7 . 

Grasping the mathematical and physical aspects of non-Markovian systems, especially with initial SS correla- 
tions, has proved to be a tough road. Nevertheless, there is a great deal of progress on deciding whether a system is 
non-Markovian in the recent years 7 " 11 . However, avoiding the initial SS correlations is not always possible in 
reality 12 " 14 . Working with initial correlations in practice has proved to be much trickier than in theory. This is 
because the presence of correlations do not allow for a clear definition of the state S independent from the state of 
S and vice versa. Physical systems are complicated and have many additional degrees of freedom that are not of 
experimental interest. Yet these extra degrees of freedom interact with the degrees of interest leading to correla- 
tions. Therefore initially uncorrelated SS state is often an approximation. 

In theory of open quantum systems, discrete quantum transformations are described by the dynamical map 
formalism 1516 : B[p s ) = pf. The dynamical map can be thought of as coming from the contraction of SS unitary 
dynamics. Let us write the state of SS as 



P 



(1) 



where yf £ is the correlations matrix 17 . The dynamical map is the mapping from the initial states of S to the final 
states of S y resulting from unitary dynamics of the SS state 



B(p s )=pf = t r£ [Up S£ tf] 

= tr f [Up s ®p £ tf] +tr £ [Ux S£ tf] 
= B CF (p s ) + B^, 



(2) 
(3) 
(4) 



where B cv is a completely positive map and B aff is the affme correction term due to the initial SS correlations. This 
means that B may not a be completely positive map when x SE # 0, nevertheless it fully describes the dynamics of 
S 19 . However, to determine such a map experimentally would require preparing different states of S while keeping 
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the SS correlations fixed. Such preparations are not operationally 
feasible because altering the state of S will also alter the SS correla- 
tions. Therefore, a nonpositive dynamical map is not an operation- 
ally meaningful quantity. 

The operational approach to quantum dynamics relies on the fact 
that quantum theory is a theory of preparations and measurements. 
The experimental method to determine a dynamical map corres- 
ponding to a quantum process is called quantum process tomography 
(QPT) 20 ' 21 . It is the central tool in determining a discrete quantum 
process; e.g. quantum gates 22 " 30 or chemical reactions 31 " 33 . To see the 
difference between QPT and dynamical maps let us review the four 
basic steps necessary to carried out QPT 34,35 : 

(i) Input states that span the space of S are prepared. 

(ii) The input states are sent through the process. 

(iii) The corresponding output states are determined by quantum 
state tomography. 

(iv) The knowledge of input states, the corresponding output 
states, and assuming linearity completely determines the 
process. 

Let us denote input states as P and output states as Q. The first step 
of QPT is state preparation. A preparation procedure takes an 
unknown state of S to a known state of S. Mathematically, it is 
described by a completely positive map acting on the system 36 . For 
instance, consider a set of preparations that project S into pure states: 

[^ (m) (x)X] (p S£ ) = P^®p^ m \ Since P {m) is a pure state, the post- 
preparation SS state is fully uncorrelated, where the con- 
ditional state of 8. X is the identity operator acting on S, as we assume 
that the preparation procedure only acts on S and not 8. We will 
discuss the implications of relaxing this assumption in Discussions. 
Lastly, if the preparation is not trace preserving, it should be divided 

tr ^4 (m) (p 5£: )J for normalisation. 

The SS evolution, after the preparation yields the output state: 

qM =tr € [u [A {m) ® j] (p se ) tf] (5) 

=tx e [upW®p e K m) rt]. (6) 

The key difference between the dynamical map in Eqs. (2) and (5) is 
the act of state preparation. Because dynamical maps of do accom- 
modate state preparation, they are not operationally defined. In the 
presence of initial SS correlations, state preparation affects the state 
of 8 in a nontrivial manner. That is, the state of 8 in Eq. (6) is 
conditioned by the choice of the preparations. 



b cp 




Figure 1 | Standard quantum process tomography. At the beginning of 
the experiment the system- environment state is uncorrelated. A 
preparation (A) is made on the system and the corresponding output state 
Q is observed. This process is described by the completely positive map of 
Eq. (4), which is a function of initial state of environment and the unitary 
dynamics. It maps the initial states of the system to output states Q. 



In deriving the standard QPT procedure it is implicitly assumed 
that the initial state of SS is uncorrelated 37 , i.e., the state of 8 is 
thought to be a constant of the problem. When that is the case, the 
state of 8 in Eq. (6) is not conditioned by the preparation procedure. 
In this case the derived map for the process is completely positive and 
is the same as the completely positive dynamical map in Eq. (4). See 
Fig. 1 for a graphical illustration. In the presence on initial SS cor- 
relation, the conditional state of 8 will be different for each prepara- 
tion, and the assumption of linearity in step (iv) of QPT is violated, 
i.e. the map is a function of the preparation procedure. Such maps are 
nonpositive, nonlinear, or simply put nonsensical 34 ' 35 . 

It then begs the question, can we determine the dynamics of a 
system that is initially correlated with £? This is an important ques- 
tion for two reasons: First, there may be physical system of interest 
that may have initial correlations. Is it possible to study their 
dynamics? Second, for foundational reasons we may care to know 
what are the limitations in describing the dynamics of physical sys- 
tems. A partial solution to these questions was given in 34 ' 38 . In this 
article we show that not only complete dynamics of initially corre- 
lated system can be determined, we can also determine the contri- 
bution due to initial correlations. 

Results 

A map on a map. QPT is performed out by noting how input states, 
that span the space of S, map to output states. The key insight in what 
follows is that it is not the input states of S that are relevant, rather it 
is the preparation procedures itself, i.e., the preparation map A. For a 
ds dimensional system there are d 2 s linearly independent states that 
span its space. However, there are d\ linearly independent operations 
(preparations) that span the space of preparations. If we determine 
the corresponding output states for a set of linearly independent 
preparations then by linearity we have can predict the output state 
for any preparation. Let us denote this map as Af-map. 

The form of Af-map arises naturally when considering the whole 
process in physical terms: At the beginning of the experiment SS is in 
an unknown (correlated) state, p se . The system is prepared into a 
known input state by the preparation procedure A, followed by a 
joint unitary dynamics. The output is given by tracing over the envir- 
onmental degrees of freedom: 

Q = tr e [U{A®l](p S£ )tf]. (7) 

We want a map acting on the preparation map A and yielding the 
output state Q: M (A) -> Q. Then Af-map is everything on the right 
hand side of Eq. (7) that is not A. The expression for Af-map in 
terms of matrix indices is 

Mrr'r"-ss's" = U r e,r'aPf>? a ,s» ^ se ,s' p' M 

Above a sum over repeated indices is implied. Af-map is a 'super 
super- operator' that acts on the super operator A. Af-map is a 
d s x d s tensor, which is contracted with a preparation A, a 
d 2 s x d 2 s tensor, yielding the output state Q, a ds x ds matrix. In term 
of matrix indices, the action is as follows: 

Qrs = Mrr'r";ss's" (Av'jsV') • ( 9 ) 

Again, a sum over repeated indices is implied. In Methods, a full 
derivation for Af-map in the last equation is given. See Fig. 2 for a 
graphical illustration of Af-map. 

Note that, in standard quantum process tomography state of 8 is a 
constant of the process, here it is the initial SS state that is the 
constant of the process, i.e., it is a fixed quantity. Physically, the 
constancy of p SE means that the experiment should be initialised 
in the same manner for every run, and then a preparation on S 
can be made. 

Al-map contains both U and p ss ; however knowing M is not 
sufficient to determine U and p S8 . As expected, it should not be 
possible to determine U and p se through measurements and 
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Figure 2 | Quantum process tomography with M -map. At the beginning 
of the experiment the system- environment state is correlated. A 
preparation is made on the system and the corresponding output state Q is 
observed. This process is described by the completely positive map M, 
which is a function of the initial system- environment state and the unitary 
dynamics. The At -map takes preparations A to output states Q. 

preparations on the system alone without access to the environment. 
Conversely, M-map contains all information necessary to fully 
determine the output state for any preparation of S. The advantage 
of dealing with the A^-map is that we have separated the preparation 
procedure from uncontrollable dynamical elements and the initial 
conditions. A^-map contains all of the dynamical information for 
the system and in the next section we will extract some of this 
information from the M-map. First let us mention some properties 
of M-map derived in Methods: Its action on a mixture of prepara- 
tions is linear, it preserves trace, it preserves Hermiticity, and it is 
completely positive. 

In Methods we show that M-map can be experimentally deter- 
mined by making a set of linearly independent preparation of the 
system. This is similar to what one has to do in standard QPT. In 
standard QPT a linearly independent set of states are fed into the 
process and the corresponding outcomes are observed. Knowing the 
inputs and the outputs the standard process map is determined. The 
difference here is that a linearly independent set of preparations are 
fed in to the process. This is of our major result of this paper: We have 
given a prescription to determine the dynamics of a system in an 
operational way, i.e., a mapping from preparations to output states. 

Quantifying initial correlations. The M-map contains the 
dynamics of the system before any preparation is made on the 
system. It is a function of the initial state p SE state as well as the 
SS unitary transformation. A^-map is a tensor, taking its trace with 
respect to the indices that belong to the initial state of S we can obtain 
the dynamics of the system as if the initial correlations we absent. 
Using this with the knowledge of the initial state of S, in Methods we 
show that from M we can derive another matrix, 



£>rr'r"]ss's" ~ ^re,r'aPr" s " Pap^se,s' p' 



(10) 



Matrix C is fully determinable from M-map and the two are the 
same when there are no initial correlations. We will call the difference 
between M and £, JC = M — £, the correlation-memory matrix: 



JC rr 'r";ss's" ~ U re f aXftf^r pU* e ^ 



(ii) 



Since M contains p SE and C contains p s ®p 8 , the difference 
between the two is a function of only x S£ - The action of the 



correlation-memory matrix on a preparation yields 

X S A {t)=K{A)=Xr £ [U[A®l\{l S£ )U% 



(12) 



which is the coherence coming into the system from the initial 
correlations. For non-Markovian dynamics the future state of 
S may depend on the initial SS correlations. This is the 



non-Markovian 'memory' due to the initial SS correlations and it 
is a key feature of non-Markovian dynamics 7 . 

The correlation-memory matrix is an important result for study- 
ing non-Markovian systems. It is an operational way of measuring 
the information that lows into S due to correlations at the time of the 
preparation. Once M-map is determined, we have the full know- 
ledge of the dynamics of S that is due to the initial correlations. The 
correlation-memory matrix provides quantitative information 
about the initial correlation and it is more than a witness for initial 
correlations 12 . 

Operational meaning of not-completely positive maps. For the 

special case, when the preparation is chosen to be the identity 
map, we get pure dynamics of the correlation-memory matrix 

JC(l)=tr B [Ux S£ U^], (13) 

which is the reduced dynamics of SS correlations. This is exactly £> aff 
in Eq. 4. From M-map we can determine matrices C and JC. In turn, 
from C we can get B cv (See Eq. (28)) and from JC we can get £> aff , and 
together they give us B of Eq. (2), which can be a not-completely 
positive map. This gives not-completely positive maps an operational 
meaning. 

Discussion 

A^-map is the result of a quantum process tomography procedure 
for initially correlated system- environment states. It is acts on the 
preparation of the initial state of the system, and only contains 
dynamical information. We study the properties of Al-map, show- 
ing it to be linear, preserving of trace and Hermiticity, and completely 
positive. Dynamical information about the evolution of the initial 
correlations can be retrieved from A^-map, in the form of the cor- 
relation-memory matrix JC. A^-map allows us to determine the 
output state for any preparation of the system, while the correla- 
tion-memory matrix JC provides a quantitative expression for the 
coherence due to the initial correlations. 

An important question is when is M-map relevant? Clearly, when 
S and S are initially uncorrelated then JC will be zero. Alternatively, 
just the presence of initial SS correlations does not warrant for M- 
map. Suppose x S£ ^0 but tYg[Ux ss U^] = 0, then the completely 

positive map of Eq. (4) would suffice to describe the dynamics cor- 
rectly for any preparation of S 18 . 

One downside to M-map is that it requires a lot of resources to 
construct. In standard quantum process tomography d 2 s input states 
are fed through the process and the corresponding output states are 
determined. To determine M map, d A s preparations are necessary, 
which is a significant growth over the standard procedure. Therefore 
an efficient way, such as compressed sensing 39,40 , to determine this 
map is desirable. This should be possible, as determining M-map is 
equivalent to carrying out ds standard quantum process tomography 
procedures. 

Another limitation that faces the procedure is the assumption that 
the preparation acts only on the system and not on the environment. 
This assumption is crucial, as we are mapping from the set of pre- 
parations on the system to the corresponding output states. If this 
assumption fails, then we would need to make a set of preparations 
that span the space of operations on the combined system- 
environment space. However, the environment can be arbitrarily 
large and we do not have any control over it. Therefore the tools 
given in this article may not be valid when the preparation affects the 
environment directly. When the preparation procedure acts on S as 
well as S, the positivity of A^-map may be affected. Note that, as long 
the effect of all preparations on S is a constant for then our prescrip- 
tion remains valid. 

Lastly, since M-map contains all dynamical information, we 
are able to construct B cv of Eq. (4) from it. Similarly, from the 
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correlation-memory matrix, we can construct £> aff of Eq. (4). 
Knowing the two we can determine B of Eq. (2), which can be a 
not- completely positive map. This gives operational meaning to not- 
completely positive dynamical map as the descriptor for the 
dynamics of the system when identity preparation is made. On the 
other hand, the non- completely positive map is not experimentally 
determinable without determining .M-map. Finally, it remains an 
open question, when JC > 0, is B not completely positive? 

Methods 

The calculations in this sections are done in terms of matrix indices as .A/f -map and 
the correlation-memory matrix JC are nontrivial tensors. We use the Einstein sum- 
mation notation, i.e., repeated indices are summed over. Bipartite state of SE is 
expressed with four indices with the Latin indices belong to S and greek indices to E. 
For instance, the state in Eq. (1) has the form pf £ s p = pf s ®P £ p + X?/s/?- ^ ma P ac ti n g 

on a density matrix is written as A rr '-ss' (pf>s>) = ^ L rr j P%^S = ^rs> where A k are the 
Sudarshan-Kraus operators (see 15 ). A* is the complex conjugation of A and A rs ^> A sr 
is the transpose; together they give Hermitian conjugation. 

A'f-map. Let us rewrite the generalised process equation, Eq. (7), in terms of matrix 
indices 



Qrs = U r€;r < 0 ,A r <r"- 1 s's" (pff^s" fi) ^se 



se;s'j3' 



(14) 



where the sum over e is the trace with respect to the environment. We are interested in 
the reduced dynamics of S as a function of the preparation procedures. Thus, we can 
pull the preparation map out of everything else and regard it all as a map acting on the 
preparation map: 

Qrs=[Ure 7 '*P%^U: e ^] (AVVS") 
= Mrr'r"-ss's" (A'r";s's")- 

In the last equation, the matrix M is defined as: 



(15) 



M r" ;ss's" = U r€j < a P$? a j, p U^j p . 



(16) 



Determining .A/f -map. Let {P (m) = |7r (m) )<7r (m) |} be a set of pure states that linearly 
span the space of S. There are d 2 s such matrices. That is, any state of S can be written 
as a linear sum of these pure states: p s = J2 nCmJ^ '• 

A preparation map acting on S is a d\ x d\ Hermitian matrix. Therefore, any 
matrix in this space can be spanned by a tensor product of the basis matrices 
|p(m) ^ p(")} 5 which is a basis in for d 2 s x d 2 s space of maps. There are d\ elements in 
the basis {P (m) (x) P (n) }. We can write action of one of these basis element on a density 
operator on S as 



A'- 



(17) 



It is crucial to note here that <7r (m) |7i (n) > 5 mn , as these vectors are eigenvectors of the 
basis elements {P (m) } that do not commute. 

These preparations are can be thought of as a projection followed by a rotation. 
Action of any map on space of S acting on the SE state can be expresses as a linear 
sum 



(18) 



= ^a<'"")| n ("))( n ('")|p S£ h ( "' ) )( ,t< " ) l ( 19 > 

mn 

= ^2 u. [mn) P {n) ®tr S [p S£ P H ] (20) 

mn 

= ^a (mn y w) P (n) (g)/l( m ), (21) 

mn 

where a (mn) are the coefficients that determine A in terms of |.4/ m "^ |. p E l(m) is the 

conditional state of the E andp( m ) = tr [P^ p s ] is the probability for the outcome P (m) . 
Knowing the output states corresponding to each of these inputs, 



Q( nui >=tr £ [UpW®p £ K m) U% 



(22) 



along with the success probabilities p (m \ for all m, n, is enough to predict the output 
state for any preparation: 



: ^ ot {mn) p {m) tr £ [l/P (M) <g>/K m ) 17*] 

mn 

^ a (»y»)Q(-). 



(24) 
(25) 



.A/f -map can be determined choosing A^ mn \ followed determining the corresponding 
Q(mn) an( jp(«) ) an( j standard inversion techniques 37 . Note that any other set of linearly 
independent preparation can be linearly mapped to the preparations given in Eq. (17), 
and therefore will suffice. 

Determining all Q (mn) is done by quantum state tomography. This is equivalent to 
carrying out ds standard QPT procedures, one each p £ \( m \ Additionally measuring 
p (m) is equivalent to doing quantum state tomography of p s . 

Before moving on a simple example may be useful. For one qubit, we may take the 
following projectors as a linearly independent basis: 

P {1) = \(I + o>), P {2) = \(I + oi), P {3) = i(I+«r 3 ), P (4) = i(I-ai). 

Note that, this is a linear but not a convex decomposition: 

p(5) = I (n + (7 2 ) =P (1) +P (4) -P (2) . The eigenvectors of P (1) , P (2) , P (3) , and P (4) 
are |* + > = -J= (|0> + 11», |>> + > = (|0> |z + > = |0>, 

and \x — ) = —=(\0) — 11)) respectively. Using these eigenvectors we can write basis 
v2 

elements for the maps that operate on the space of one qubit. For instance, 

A {hl) (p S£ ) = \x^)(x^\p S£ \x^)(x^\=p^\x^)(x^\®p^ x+ \ 

A {XA) (p S£ ) = \x-)(z+\p S£ \z + )(x-\=p^\x-){x-\®p £ ^ + \ 

^(4,2) (p^) = | r+ >(^_|p^|^_>( y - h | = p(--) b - h > 0/ - h | ( x)^l(--), 

and so on. 

Detecting initial correlations. The initial state of the system is labeled by indices f 
and s". Tracing over everything else we can find the initial state of S (before 
preparation) from A4-map: 

1 s 

— d rs d r > s >Mrr>r";ss>s" = P^^ ■ (26) 

d s 

This is, of course, attainable by doing state tomography at the beginning of the 
experiment, by measuring the values of p (m) from last section. 
Next, let us the trace over the system indices r" and s" 

8r*fMr,r"tffi = U^ a pZ p U^p = &%.j. (27) 

The last equation is exactly the dynamical map in the absence of initial correlations, 
given in Eq. (4). In other words, in the absence of initial correlations, QPT would yield 
this map. 

This means, even though the M -map contains the information about uncorrelated 
SE state and the correlations separately. Consider the following matrix composed of 
the matrices in Eqs. (26) and (27) 

r — n CP n s 

*~>rr'r";ss's" ~ ° rr' :ss' P r" s" 

(28) 

= Ure,r'aP%'P £ a pU^ p . 

The last equation is similar to the expression for the A4-map, except the state of the 
system and the state of the environment are uncorrelated. 
Writing the state of SE in M-map in terms of Eq. (1), we get 



Now we can define the correlation-memory matrix as 

/C rr V' ;ss ' s " =Mrr l r";ss's" —■Crr'ru-^s" 
= U rej i a Xf £ a - S " p U* e j p . 



(29) 

(30) 
(31) 



Q = tv s [UA(p s£ )uT] 



(23) 



Properties of M. Linearity. Mathematically, A4-map acts on the preparation map 
just as the dynamical map acts on a density operator. In fact, we are not varying the 
initial state of the system, rather the preparation procedure on that state. Therefore 
the linearity of quantum mechanics is preserved for the A4-map acting on different 
preparation procedures, i.e. 

M[ ai ^ (1) +a 2 ^ (2) ] =a l MA [l) +d 2 MA {2) . (32) 

This is very much like the dynamical maps action on mixtures of states. Furthermore, 
if we show that the .A/f -map preserves trace, Hermiticity, and positivity on its domain 
then all of these properties will be preserved on the state space. In other words for any 
preparation, A^ that preserves trace, Hermiticity, and positivity, the action of the 
M-map on it will yield an output state, Q (m) , that is unit-trace, Hermitian and 
positive. 
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Trace preservation. Let us start with the trace of M with respect to the final indices r 
with s: 

tr {rs) [M]=5 rs Mrr'r";ss's" 

Since tfU=Zre U^ p U r e,r>* = I> then 

tr (rs) [M] = d^d^pf^ ip = mp S . (34) 

A preparation acting on the above matrix will yield 

tr irs) [M](A)=tr[A(p s )]=l. (35) 

The implication being M preserves the trace of A(p s ) . As long as the preparation is 
trace a preserving operation we get a unit-trace matrix for the output state. 

Hermiticity preservation. As with the case of general quantum operations, matrix Ai 
is Hermitian. This is easy to see by taking the complex conjugate of matrix M, 

*pe ( 36 ) 

= M ss >s";rr"r<- 

The complex conjugate of M is not only the transpose of M, but each element of M 
is also transposed. Hence Ai is a Hermitian matrix. 

Positivity of Ai-map. The A4-map is composed of a unitary matrix operating on a 
density matrix. Then we can take the square root of the density matrix to get 

(M r r>r";ss>s») = I/*/ a \/p^ r »a,ay ^/l^^ 

* (37) 

= M r/r» M s^:s, 

where M = Uy/ p S£ and ji = aye. We have written the -map in operator sum 
representation, hence it is completely positive. Where M At are the Sudarshan-Kraus 
operators 1516 . This means, the .A/f -map acting on any preparation procedure will lead 
to a physical state. This was not the case when a standard QPT procedure is carried 
out on initially correlated SE states. The action of A4-map can now be written as 

/< 

The properties shown above are precisely the conditions for a generic quantum 
operation to preserve trace, Hermiticity and positivity. Therefore A4-map preserves 
the attributes on the preparations, which in return will preserve these attributes on the 
states. 
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